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On the relation between virial coefficients and the close-packing of hard disks 
and hard spheres 
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The question of whether the known virial coefBcients are enough to determine the packing fraction rj,^ at 
which the fluid equation of state of a hard-sphere fluid diverges is addressed. It is found that the information 
derived from the direct Pade approximants to the compressibility factor constructed with the virial coefficients 
is inconclusive. An alternative approach is proposed which makes use of the same virial coefficients and of the 
equation of state in a form where the packing fraction is explicitly given as a function of the pressure. The 
results of this approach both for hard-disk and hard-sphere fluids, which can straightforwardly accommodate 
higher virial coefficients when available, lends support to the conjecture that 7700 is equal to the maximum 
packing fraction corresponding to an ordered crystalline structure. 
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I. INTRODUCTION 

The virial expansion of the equation of state of a fluid 
is an expansion in powers of (usually) the number density 
p that reads 



Zip,T) = l 



J=2 



(1.1) 



where Bj are the virial coefficients and Z = p/phsT 
is the compressibility factor, with p, ks, and T being 
the pressure, the Boltzmann constant, and the absolute 
temperature, respectively. This expansion was originally 
introduced by Thiesen in 1885 (Ref. [l| in connection with 
the equation of state of fluids at low densities. Soon 
afterwards, and apparently independently, it was used 
by Kamerlingh Onnes^ (who named the coefficients in 
the expansion as virial coefficients), in order to provide 
a mathematical representation of experimental results. 
The virial series was later proven to arise naturally in 
rigorous derivations in statistical mechanics^ whereby the 
Bj turn out to be related to intcrmolecular interactions 
and are in general functions of temperature. 

In the case of hard-core systems such as hard disks or 
hard spheres)^ the virial coefficients do not depend on 
temperature. In particular, the value of the second virial 
coefficient for hard spheres of diameter a in d dimensions 
is B2 = where Vd = in/4f/'^/T{l + d/2) is the 

volume of a d-dimensional hard sphere of unit diameter, 
a result first derived for three-dimensional hard spheres 
(d = 3) by van der Waalsi^ Analytical expressions for 
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TABLE I. Known virial coefficients^ for a hard-disk fluid 
{d — 2) and a hard-sphere fluid {d — 3). 
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bj {d = 2) 


bj {d = 3) 
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2 
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3.12801775... 


10 
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4.25785445 . . . 


18.36476838 . . . 


5 


5.3368966(2) 


28.22451(26) 


6 


6.36296(13) 


39.81515(93) 


7 


7.35186(28) 


53.3444(37) 


8 


8.31910(44) 


68.538(18) 


9 


9.27215(90) 


85.813(85) 


10 


10.2163(41) 


105.78(39) 



B3 and B4 are also available in the literature^"— but 
higher virial coefficients must be computed numerically 
and, since this represents a non trivial task, up to now 
only values up to the tenth virial coefficient have been 
reportedFii"— 

The virial expansion for d-dimensional hard-sphere 
systems is often expressed in terms of the packing frac- 
tion 77 defined as 



77 = Vdpa 



(1.2) 



Hence, for these systems the virial expansion of the com- 
pressibility factor is given by 



ziv) = 1 



i=2 



(1.3) 



where the (reduced) virial coefficients bj = Bj/ (wdcr 
are pure numbers. For one-dimensional hard rods, hj = 1 
and so Z{ri) = 1/(1 — 77). The presently known values^ 
of the virial coefficients coefficients for hard disks [d = 2) 
and hard spheres {d — 3) are given in Table HI 

There are a number of controversial and open issues 
related to the virial expansion of hard-disk and hard- 
sphere fluidsi^ To begin with, even in the case that many 
more virial coefficients for these systems were known, the 
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truncated virial series for the corresponding compress- 
ibility factors would not be useful in principle for pack- 
ing fractions higher than the one corresponding to the 
radius of convergence ?7conv of the whole series. Such a 
radius of convergence is determined by the modulus of 
the singularity of Z{r]) closest to the origin in the com- 
plex plane. While its actual value is not known, lower 
bounds are available^^i^S and existing evidence suggests 
that it derives from a singularity located on the real neg- 
ative axisi ^^'^° In fact, one of the major reasons for trying 
to evaluate higher order virial coefficients is precisely the 
determination of the value of Tyconv and of the nature of 
the singularity giving rise to it. One should add that, al- 
though the virial series diverges for 77 > yyconv, given the 
fact that r^conv seems to be located outside the positive 
real axis, one would expect that the compressibility fac- 
tor would still be well defined for 77 > 7/conv, at least in a 
certain range. 

Coming back to the controversial and open issues, even 
the character of the virial expansion (either alternating 
or not) is still unknown. So far, all the available virial co- 
efficients for these systems are positive, but results from 
higher dimensions suggest that this feature might not 
be true for higher virial coefficients 1^ Finally, the evi- 
dence coming from approximate equations of state ob- 
tained through the knowledge of the limited number of 
virial coefficients via various series acceleration methods, 
such as Padc or Levin approximants, indicates that the 
freezing transition observed in computer simulations does 
not show up as a singularity in these equations of state 4i 
As a matter of fact, while it is quite plausible that Z{'q) 
presents a singularity at the freezing density rjf^— the 
virial coefficients (or even their asymptotic behavior) do 
not seem to yield any information concerning the freezing 
transition at 7^/4^ This might be related to the fact that 
lim^^^- Z{r]) = finite. 

Since the compressibility factor of hard-disk and hard- 
sphere fluids both for the stable and metastable fluid 
phases is a monotonically increasing function of the pack- 
ing fraction;^ one may reasonably wonder at which pack- 
ing fraction 77 = rjao the analytical continuation of the 
compressibility factor diverges, namely one wants to find 
the value of 7700 such that 

lim Z{'q) = 00. (1.4) 

Clearly, rjoo may not be bigger than the maximum pack- 
ing fraction T^max that is geometrically possible (the so- 
called Kepler's problem). Therefore, one must have 

Vconv !i Voo '7max- (1-5) 

The value of T^max, at least for not too high dimension- 
alities d, corresponds to an ordered crystalline structure. 
In Table |TT] we provide the values of 77max for d = 2-8. 

More than a decade ago, by studying the singularities 
of Padc approximants constructed from the virial series 
for hard disks and hard spheres, Sanchca^ came to the 



TABLE 11. Values of 77n,ax for d = 2-8.- 



d Exact Numerical 

2 0.9069 

3 lnV2 0.7405 

4 ^Tv^ 0.6168 

5 io'^^^ 0.4652 

6 T34^^V3 0.3729 

7 j^tt'^ 0.2953 

8 :^7r* 0.2537 



conclusion that such singularities were related to crys- 
talline close-packing in these systems. Other authors^"— 
have also conjectured that 

Voo = ?7max- (1-6) 



It is interesting to note that this conjecture was already 
suggested by Korteweg^ and Boltzmann^ in the late 
1800s. 

On the other hand, the conjecture (|1.6p has not been 
free from criticism^ and some authors^"— have con- 
jectured that rjao = r/rcp, where rjy-cp (approximately 
equal to 0.82 and 0.64 for hard disks and hard spheres, 
respectively^) is the random close-packing fraction. For 
a thorough account of proposed equations of state, in- 
cluding those enforcing tJoo = ?7max or 7700 = ?7rcp, the 
reader is referred to Rcf. |70i. 

In his original study, Sanchez used a Fade analysis 
with the then most recent available values of the first 
eight virial coefficients given by Janse van Rensburgj^ 
More recently. Sanchez and Lee^ found that, using the 
ninth and tenth virial coefficients^ for a hard-sphere 
fluid, the corresponding Fade approximants remained fi- 
nite at 77 = r/inax- However, taking into account such 
coefficients, they constructed a new approximate equa- 
tion of state that diverges at that packing fraction. In 
this regard, the question arises as to whether using the 
available information on virial coefficients one can reach 
a more definite conclusion. More generally, one could ask 
whether there may be a systematic method to improve 
the estimation of 7700 as more virial coefficients become 
available. 

The paper is organized as follows. In Sec.|lT]we perform 
an analysis of the singularities of the compressibility fac- 
tors obtained from the different Fade approximants con- 
structed from the known virial coefficients of the hard- 
disk and hard-sphere fluids. This is followed in Sec. Hill 
by the introduction of inverse representations of the virial 
series and their connection with the computation of rjao , 
the results being presented in Sec. lIVl The paper is closed 
in Sec. |V] with some concluding remarks. 
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II. PADE APPROXIMANTS AND SINGULARITIES IN 
THE COMPRESSIBILITY FACTOR 

In general, given the series 



(2.1) 



the notion behind a Padc approximant is the replacement 
of the series by a ratio of two polynomials, namely 



(2.2) 



where, without loss of generality, one may take /3o = 1- 
The remaining N + M + 1 coefficients are chosen in such 
a way that the Taylor series expansion of P^^ (z) exactly 
yields the first N + M +1 terms of the power series S{z). 
On the other hand, the difference iV — M of the degrees 
of the two polynomials may be taken at will. A more 
detailed treatment of Padc approximants may be found 
in Ref. 

In the usual application of Fade approximants for the 
equation of state of hard-sphere fluids, the compressibil- 
ity factor Z is approximated by 



ziv) 



pM 



iv) 



(2.3) 



where, since lim^_>.o Z{t]) = 1, one sets ao ~ 1. The 
other N + M coefficients {aj,j = 1, . . . , M} and {(3j,j = 
1, . . . , N} are determined with the values of the virial co- 
efficients {bj,j ^ 2, . . . N + M + I}. For a given value of 
g = N + M, there are in principle g independent approxi- 
mants (corresponding to = 1, . . . , <?), because the case 
{N,M) = {0,g) is not a Fade approximant but rather a 
truncated series. We will use the term 'Fade degree' for 
g = N + M . Given the known values of bj (see Table U), 
the highest value of the Fade degree is 5 = 9. The Fade 
approximants (|2.3p are not only used to get the equa- 
tion of state but also to predict values of unknown virial 
coefficients 

Both the radius of convergence r^conv of the virial series 
and the value of rjoo when Z is approximated by (|2.3p 
are determined by the nature of the singularities of the 
approximant. The modulus of the complex root of 1 -I- 

Ej!=i l^jV'' closest to the origin gives r^conv: while r^oo is 
identified with the smallest real positive root. 

Due to our interest in obtaining the value of rjoo, we 
have carried out a systematic analysis of the location of 
the real and positive pole of {i]) closest to the ori- 
gin. For each Fade degree g < 9 we have considered 
all the associated g Fade approximants, namely a total 
of Eg=i = 45 approximants, and located their poles. 
In some instances, the corresponding approximant either 
has no real and positive poles or these are located be- 
yond r/ = 1. In the rest of the cases, if more than one 
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FIG. 1. Value of the real and positive pole closest to the 
origin for the different Fade approximants of the form p. 31) as 
a function of the Fade degree g = N + M. The number placed 
next to each circle indicates the degree A'' of the polynomial 
in the denominator of the corresponding Fade approximant. 
The horizontal lines correspond to 77 = ?7max. 



real and positive pole existed in the range < 77 < 1, we 
have focused on the smallest one of them. Moreover, we 
have discarded a pole when it practically coincided with 
a zero of the approximant. The results of this analysis 
are shown in Fig. [TJ 

As can be seen in this figure, only 29 {d = 2) or 23 
(d = 3) out of the 45 approximants possess a real and 
positive pole in the range < 77 < 1. The Fade approx- 
imants with A = 1 have a real pole at ?7poie = bg/bg+i 
and this is the only pole for a few values of g. Also 
worth pointing out is the fact that the real and positive 
singularity closest to the origin exhibits an important dis- 
persion for a given Fade degree g, something particularly 
noticeable for <? = 8 and d = 3. Moreover, for <? > 5 the 
value of the poles may even be higher than 77max, which 
is clearly unphysical. 

In order to check on the robustness of the method of 
analysis, we have carried out the same analysis for the 
cases of (7 = 8 and 3 = 9, this time varying the values of 
the two last known virial coefficients by 5%. More specif- 



4 



1.00 

0.95 

0.90 

0.85 

0.80 

1.0 
0.9 - 
0.8 - 



0.7 
0.6 ^ 
0.5 
0.4 



3& 



1* 

29. 

70 



1|6 



2o 
70 



5. 

It 6. 
7" 4 

2q7 



50 



H 1 1 ^ 



d=2 



d=3 



10 



III. INVERSE REPRESENTATIONS OF THE VIRIAL 
SERIES 

We have seen that the use of the direct Pade approx- 
imants of the comprcssibihty factor Z is not reUablc for 
the purpose of determining 7700 ■ This is perhaps not sur- 
prising since, as McCojs^ has pointed out in his recent 
book (and we quote), "there is no reason to expect that 
simple model equations of state which use the density 
as the independent variable will capture the true physics 
of hard particle systems (or indeed of any real system)." 
Therefore, an alternative approach is called for. We be- 
gin by introducing an equivalent form of the virial series 
p.3p . namely 



P 



0=2 



(3.1) 



Note that p = -qZ . Following an idea of Sanchez,— we 
may formally invert the series in ()1.3p . leading to 



(3.2) 



,?(Z) = ^e,(Z-l) 



Here we propose a similar inversion of the series (|3.ip in 
the form 



">]{¥) = P ^ 



E 



c.p' 



(3.3) 



FIG. 2. Value of the real and positive pole closest to the origin 
for g = 8 and g = 9. Filled circles correspond to the values ob- 
tained from the use of the coefBcients 69 and &10 given in Table 
in while the open circles have been obtained using the values 
(69, hw) = (8.80, 9.71) for d = 2 and (69, 610) = (81.52, 100.5) 
for d = 3. The number placed next to each circle (in ital- 
ics for open circles) indicates the degree A'' of the polynomial 
in the denominator of the corresponding Pade approximant. 
The horizontal lines correspond to r^max- 



ically, instead of the values of 69 and 610 given in Table 
m we have taken (&9,&io) = (8.80,9.71) for d = 2 and 
(&9,&io) = (81.52, 100.5) for d = 3. The consequence of 
such replacement is an important variation in the charac- 
ter and location of the poles, as shown in Fig. [21 includ- 
ing the appearance of new poles and the disappearance 
of earlier ones. An exception is the pole r/poic ~ &9/&10 
corresponding to = 1 and g = 9, which obviously is 
not affected by a common factor multiplying both 69 and 
610. Given this generally great sensitivity of the poles on 
the values of the virial coefficients, it does not seem wise 
to rely on the previous analysis in order to get good esti- 
mates of ?7oo. Instead, in Sec. lIIII we will follow a different 
route. 



The coefficients {cj} and {cj} arc (nonlinear) combina- 
tions of the virial coefficients {bj}. In fact, the determi- 
nation of Cj and Cj requires that one previously knows 
j and j -f 1 virial coefficients, respectively. Notice that, 
while going from Eq. (|1.3p to Eq. p.ip or vice versa is 
a trivial task, the same is not true when considering the 
inverse developments (|3.2p and p.3p . Further, the con- 
dition p.4p to obtain ?/oo is equivalent to 



lim ij{Z) = lim ri{p). 

p — ^00 



(3.4) 



The form p.2p was the one that Sanchez used^ to esti- 
mate the value of 7700. However, the true thermodynamic 
variable, along with the density (hereby represented by 
77), is the pressure (hereby represented by p), and not 
the ratio pressure over density (given by Z). Therefore, 
in our view the form p.3p has a clearer physical mean- 
ing than p.2p . In this respect, it is interesting to note 
that Hamad^iS proposed for hard-disk and hard-sphere 
fluids an approximate equation of state of the form 



vip) 



P 



l + b2p- 



(3.5) 



with ki — (62 — by,)/{k2 — fcs) and where ^2 and fcs were 
obtained by a fitting method. According to Eq. p.Sp . 
^00 = [&2 ~ ki \a{k2 / k^)]~^ , which turns out to be close 
to the respective random close-packing valuesi^ 
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Since bg and 610 are of the order of 10, it is clear that nei- 
ther cg nor cio are strongly affected by the precise values 
of 6g and 610. This effect is even much more pronounced 
in the case of hard spheres, where 



eg = 2.78 X 10^ (1 - 3.60 x lO-'^&g) 



(3.10) 



cio 



-2.15 X 10^ (1 - 2.05 X 10^^69 + 4.66 x 10~^&io) • 

(3.11) 

Taking into account that both bg and 610 are of order 
10^, one sees that the actual values of cg and cio are very 
weakly influenced by bg and 610. The same process can 
be repeated for cn and C12, this time using the values of 
bg and &10 listed in Table U but leaving 611 and &12 free. 
The results for hard disks arc 



cii = 2.72 x 10^ (1 - 3.67 x 10"^6ii 



(3.12) 



FIG. 3. Semi-logarithmic plot of the absolute value of the 
inverse virial coefficients Cj for j < 10. The dashed lines 
represent the exponential fit \cj\ ~ C9{\cio\/cgy~^ . 



There is an additional argument in favor of expressing 
the equation of state in the form rj ~ fjip)- For one- 
dimensional systems of particles interacting through an 
arbitrary potential ip{r) (restricted to nearest neighbors) 
there is an exact statistical-mechanical solution.—"— In 
this solution, the equation of state is given by2^ 



no{p/kT,T) 
fliip/kT, T)' 



where 



dr r^e 



Sf,-xr-ip(r)/kT 



(3.6) 



(3.7) 



Thus, at a given temperature, the density is an explicit 
function of the pressure. In general (the hard-rod case is 
an exception) , it is not possible to invert Eq. (|3.6p to ex- 
press the pressure as an explicit function of density or the 
density as an explicit function of Z . The functional form 
p.6p may be easily extended to the case of mixturesi^ 

We have obtained the values of the inverse virial coef- 
ficients Cj for 2 < j < 10 and Cj for 1 < < 9 from the 
knowledge of the first ten virial coefficients listed in Ta- 
ble HI Here we give some details about the coefficients Cj. 
It is observed that they alternate sign (c2 ~ —62 being 
negative) and their absolute values grow almost expo- 
nentially with j . The results are graphically displayed in 

Fig.m 

For further use, it is interesting to provide the expres- 
sions of the coefficients cg and cio when 52-&8 are given 
by Table |T] but bg and 610 remain free. For hard disks, 
one gets 



Cg = 2.86 x 10^ (1 - 3.49 x lO'^^^g 



(3.8) 



cio 



-8.97 X 10^ (1 - 2.45 x 10~^69 + 1.11 x 10~'*6i 



C12 = -8.53 X lO'' (1 - 3.05 X + 1.17 x lG"^6i2) . 

(3.13) 

Analogously, for hard spheres one finds 

cii = 1.68 X 10® (1 - 5.96 X lO^^fen) , (3.14) 



C12 -1.32 X 10^ (1 - 3.92 X 10"**6ii + 7.55 x 10-^%V2) 

(3.15) 

Therefore, as j increases, the values of Cj become less and 
less sensitive to the actual values of the virial coefficients 
bj-i and bj, this effect being much more important for 
hard spheres than for hard disks. This implies that one 
could get good estimates of cn and C12 (especially for 
hard spheres) even with poor estimates of the unknown 
virial coefficients 611 and 612- We will come back to this 
point later. 

Once we have introduced the representations (|3.2p and 
p.3p . the next step is to compute their corresponding 
Fade approximants. These read 



{Z-l)P^- 



and 



\Z-l) 



(3.16) 



(3.17) 



Now one can determine 7700 using Eq. p.4p . Since rjoo 
must be finite and different from zero, we must restrict 
ourselves to diagonal Fade approximants, in contrast 
with what occurs with the direct approximants 
Therefore, 



r]oo{N) 



CtN-l 



N 



(3.18) 



(3.9) 



where, following the notation of Eq. (|2.2p . a^v-i and 
are the coefficients of the terms of the highest degree 
in the numerator and denominator of the approximant, 
respectively. In Eq. (|3.18p the notation T]ao{N) indicates 
that, in principle, the value so obtained may depend on 
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TV. One would expect that the true value of rjoo would 
be 77oo = limAT^oo Voo{N). 

Smce 61 = 1, m the case of Eq. p.l6p we have to 
determine 2N coefScients usmg ei,...e2N (or, equiva- 
lently, 62,. . . f'27V+i)- On the other hand, in order to use 
Eq. (|3.17p . 2N — 1 coefficients must be determined from 
C2,. . . C2N (or, alternatively, from 62,. . . fe2w). We recall 
that only ten virial coefficients arc presently known, so 
the highest N that in principle may be taken in connec- 
tion with Eqs. ()3.16|) and p.l7|) is 4 and 5, respectively. 
This restriction may be removed at the expense of includ- 
ing estimates of the first few virial coefficients beyond the 
tenth. 



IV. RESULTS 

The results for r7oo(-^)i derived from the different Pade 
approximants as given by Eq. (|3.16p (with TV = 1, 2, 3, 4) 
and the use of Eq. (|3.18p arc shown in Fig. |4l Those 
pertaining to Eq. ((XT7)) (this time with TV = 1, 2, 3, 4, 5) 
appear in Fig. [5l 

As clearly seen from these two figures, the behavior 
of the inverse representations is much more regular than 
the one displayed in Fig. [T] This strongly supports the 
notion that they are more adequate if one wants to get 
a reliable value for 7700- Note that the results derived 
from Eq. (|3.17p . (c/. Fig. [5]) show a smooth behavior and 
that the estimate rjooiN) increases slowly with TV. In the 
three-dimensional case, it is not clear whether r]oo{N) 
converges to the close-packing value Jjmax for TV — > c« 

0.64. How- 



rep 



or to the random close-packing value rj, 
ever, the very good fit obtained with a law of the form 
^700 (TV) = ^ymax — ae^^^ suggests that the values ob- 
tained up to TV = 5 are not incompatible with the result 
limjv-j.00 Voci^) = 'ymax- In the two-dimensional case the 
approach of 7700 (TV) to jymax is much clearer than in the 
three-dimensional case, the value of 7700 (5) being even 



slightly larger than rj^cp — 0.82. This gives further sup- 
port to the correctness of the conjecture (|1.6p . 

Now we turn to Fig. Although the representation 
p.l6p leads to one point less in the graph than Eq. p.l7p 
and does not make use of 610, ??oo(4) is remarkably close 
to r?max, especially in the case d = 2. This could be 
just a coincidence and it is conceivable that, once bn 
became available and the point corresponding to TV = 5 
were included, 7?oo(5) would turn out to be greater than 

In fact, the fit to a law r?oo(TV) 
is much poorer in this instance. Further, we have found 
that the results for ??oo(TV) obtained from Eq. (|3.16p are 
more sensitive to variations of the values of the virial 
coefficients than those obtained from Eq. p.l7p . 

In view of the above, and besides the physical reasons 
we alluded to before, we find that the pressure represen- 
tation seems to be the most reliable one for estimating 
the true value of r]oo- On the other hand, it seems also 
clear that, strictly speaking, the knowledge of the first 
ten virial coefficients is not enough to decide whether 
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in the three-dimensional case. 



In order to look into the performance of the represen- 
tation (|3.17|) in more detail, we have carried out a test of 
robustness similar to the one made in connection with the 
poles of the conventional approximants in Fig. [31 In this 
case, however, rather than considering for each TV only 
an error of ±5% in the values of the highest virial coeffi- 
cients, we take an error of ±100% in 62 w- The results of 
this procedure are also shown in Fig. [SJ where it is clear 
that as one increases TV the effect of the error becomes 
less pronounced, to the point that for TV = 5 and d = 3 it 
becomes practically unnoticeable. This means that one 
cannot rely on the representation (j3.17p for getting esti- 
mates of unknown virial coefficients. On the other hand, 
its introduction was not made for that purpose but rather 
for obtaining good estimates of rjoo- 

In this regard, albeit on a more speculative basis, for 
the sake of going beyond the limit imposed by the value 
TV = 5 one can incorporate into the procedure estimates 
of unknown higher virial coefficients fen, 612, .... Given 
the robustness of the results as discussed above, one 
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FIG. 5. Estimates of r}oc derived from Eq. (|3.17p . The hor- 
izontal hnes correspond to ?7max- The curves correspond to 



the fit r]ac,{N) = rjn 



0.65e- 



(d = 2) and r]oo{N) 



?7max — 0.63e~ ■ (d = 3). The circles, squares, and dia- 
monds correspond to taking for &2Ar its correct value given in 
Table [H zero, and twice the correct value, respectively. 



would expect that the effect of errors in such estimates on 
VooiN) for N somewhat greater than 5 would be weak. 
The idea that this is indeed the case goes back to Eqs. 
(|3.8p - (|3.15|) . Let us specialize here to d = 3. From Eqs. 
((3T0)) and ((3lT|) one obtains 



a4 



120.505 



1 - 1.30 X 10-359 
-1.068 X 10"^6q 4 



2.64 X 10~^6g 



4.40 X IQ-'^b 



211.567 



1 - 1.30 X 10-36q 



10; 



2.65 X 10-^69 



-1.31 X 10~^5^ +4.79 X 10-% 



10 



(4.1) 



(4.2) 



for the coefficients of highest degree in the numerator 
and denominator of the Fade approximant p.l7p with 
N = 5. Since 69 ~ &10 ^ 10^, it turns out that a4 and 
/Js are hardly affected by the precise values of and 
&10. But, because of a partial cancelation of terms, this 



TABLE III. Estimated virial coefRcientsS^ fen- 
disk fluid (d — 2) and a hard-sphere fluid {d = 



■6 16 for a hard- 
3). 



bj {d ^ 2) 



{d ^ 3) 



11 
12 
13 
14 
15 
16 



11.15 
12.08 
13.03 
13.93 
14.91 
15.86 



128 
153 
182 
215 
247 
279 



influence is still weaker in the case of the ratio 0,4^ / (3^^ 
namely 

?7oo(5) ~ 0.5696 (1 -f 4.6 x + 2.5 x 10"^6| 

-3.9 X 10^^610) . (4.3) 

When replacing the known values 69 — 85.813 and 5io = 
105.78 one gets 7700 (5) = 0.5696 x (1 - 0.002), which only 
differs 0.2% from the value obtained by setting 69 = 5io 
0. 

Let us repeat the same process in the case iV = 6. It 
can be checked that Eqs. (|3.14p and p.l5|) yield 

- 5.98x10-%, 



as 



1 - 2.95 X 10-4foii" 
-6.03 X + 1.11 X 10-'^6i2) 



(4.4) 



482.486 



;;i - 5.99 X 10-'*6ii 



1 - 2.95 X IQ-^ii 
-6.98 X + 1.18 X 10-%2) , 



(4.5) 



7700(6) ~ 0.6045 (1 + 5.1 X 10"%i + 9.4 x IQ^^bl^ 



-7.1 X lO^^foi 



(4.6) 



Inserting the estimated values 611 ~ 128, 612 — 153 (see 
TableHni) we obtain 77^0(6) 0.6045 x (1-0.0008), which 
deviates less than 0.1% from the value corresponding to 
611 = 612 = 0. Considering the above as an illustra- 
tive example, it is not unreasonable to expect that, as 
increases, the value of rjao{N) becomes less and less 
sensitive to the actual values of the virial coefficients bj 
with j < 27V. 

With the previous assertion in mind, here we will go 
up to A'' = 8 taking the values of 611-616 as estimated 
by Clisby and McCoy^ from Fade approximants. These 
estimates are given in Table IIIII The corresponding re- 
sults for rjoa{N) are shown in Fig.[6l As was done in Fig. 
[51 in order to illustrate representative outcomes, we have 
also included in Fig. |6] for iV = 6, 7, and 8 the points 
obtained when the values of 611 and 612 from Table IIIII 
for d = 2 (d = 3) have an error of ±10% (±20%), those 
of 613 and 614 an error of ±20% (±40%), and those of 
615 and 616 an error of ±30% (±60%). Once more, irre- 
spective of the errors, one can not distinguish the three 
points obtained in this way with the same A'^ > 6 in the 
case d = 3. Apart from this, perhaps the most interesting 
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FIG. 6. Estimates of 7700 derived from Eq. p.l7[) . The results 
up to TV = 5 (filled circles) have been derived from the known 
values of the virial coefficients as given in Table [H while those 
for TV = 6, 7, and 8 (open circles) have been obtained from 
the estimates of the higher virial coefficients given in Table 
mil The open squares and diamonds represent the values of 
??oo(TV) when the estimates bn-bie of Table Hill are assumed 
to have an increasing error between ±10% and ±30% {d — 2) 
or between ±20% and ±60% (d = 3), as explained in the 
text. The horizontal lines correspond to 77max. The curves 
correspond to the same fits as in Fig. [S] 



feature that one can observe from Fig. [HI is that the new 
additions appear to continue the smooth trend obtained 
with the previous five points. In fact, they fall on top 
of the curve that served to fit these previous points in 
Fig. [S] Further, the last point of the curve correspond- 
ing to d = 3 is above the random close-packing value 
r^i-cp ~ 0.64. Hence, all this evidence strongly supports 
the conjecture (|1.6p . However, a word of caution is called 
for since, irrespective of the tests that we have carried out 
on the sensitivity of 7700 (TV) to errors in the values of the 
virial coefficients, the possibility that some of the esti- 
mates that we have used for the higher virial coefficients 
may be quite deficient does exist. One must recall that in 
higher dimensions there are negative virial coefficients^ 
and, although it is unlikely that for hard spheres some 



of the coefficients bu to biQ may be negative, one cannot 
rule it out with certainty. 



Since the results of this section face the limitation of 
having been derived with the (relatively) small number of 
available virial coefficients, one may reasonably wonder 
whether our proposal would stand against a more strin- 
gent test. For instance, one could examine the model of 
hard hexagons on a lattice, for which an exact solution 
exists and where the first twenty five virial coefficients 
are given by Joyce or the model of hard squares on a 
lattice, where the first forty two virial coefficients can be 
obtained from the work by Baxter et al^ While the anal- 
ysis (not shown but available upon request) carried out 
for these systems using the inverse Fade representation 
of the pressure similar to Eq. (|3.17l) is not incompatible 
with the validity of the conjecture (|1.6p . it presents large 
oscillations and is flawed by unphysical results arising 
from some defective Fade approximants. This may be 
due to different causes which we will discuss in connec- 
tion with the hard hexagon model on a lattice. On the 
one hand, continuous and two-dimensional discrete mod- 
els show important differences. For instance, in the hard 
hexagon model on a lattice there is no analytic contin- 
uation of the pressure from the low to the high density 
regime. In fact, the fluid branch and the solid branch 
both end and meet with a zero slope at the critical den- 
sity in what appears to be a continuous phase transition 
from the disordered to the ordered state. In this case, 
no metastable states are possible. In contrast, computer 
simulations have revealed both for the hard-disk and the 
hard-sphere fluids that, in the thermodynamic limit, the 
low and high density branches of the equation of state 
do not smoothly join one another. In the case of hard 
disks, it is not clear yet whether the fluid-solid transi- 
tion is first order or mediated by an intermediate hexatic 
phase. In the case of hard spheres, there exists com- 
pelling evidence that these systems present a first order 
phase transition where a two-phase region appears with 
a constant pressure as the density varies from the den- 
sity of the pure fluid p/ to the density of the pure solid 
Ps- The flat portion of the isotherm which connects pf 
and Ps is the tie line. The fluid branch for p > pf is 
metastable but can be observed in computer simulations. 
Thus, while in the case of the hard hexagon model on a 
lattice one does know that the pressure diverges at the 
close-packing density in the solid branch, the question 
of where the analytic continuation of the fluid branch 
may diverge loses all meaning. On the other hand, as 
mentioned in Sec. U for the hard-disk and hard-sphere 
fluids the existence of metastable states past the freezing 
point legitimately allows one to ask at which density this 
metastable fluid branch, when analytically continued, di- 
verges. The present results suggest that knowledge of the 
virial coefficients may indeed indicate when this will oc- 
cur for (continuous) hard-core fluids. 
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V. CONCLUDING REMARKS 

On the basis of our results, the foUowing conclusions 
and remarks can be made. Two issues have been ad- 
dressed in this paper. The first one concerns the ques- 
tion of whether the known virial coefficients (presently 
ten both for hard disks and hard spheres) are enough to 
tell us about the packing fraction r]aa at which the fluid 
equation of state, continued and extrapolated beyond the 
fluid-solid transition, has a divergence to infinity. In con- 
nection with this issue, the determination of 7700 from the 
direct Pade approximants [c/. Eq. ()2.3|) ] seems to be not 
very reliable. Given the fact that as a thermodynamic 
variable the pressure has a clear physical meaning not 
shared by the compressibility factor, the inverse repre- 
sentation (|3.17|) appears then as the natural candidate 
for such a determination, at least when restricting to 
the presently available information, namely only the ten 
known virial coefficients. 

The second issue is related to the validity of the con- 
jecture p.6p . In the case of hard disks, all the results 
support it. On the other hand, strictly speaking, for 
the hard-sphere system our results allow us neither to 
validate nor to discard this conjecture. However, the ro- 
bustness analysis carried out for hard spheres strongly 
suggests that the conjecture is true. Nevertheless, full 
confirmation must await the availability of higher virial 
coefficients. In any case, we are persuaded of the useful- 
ness of the inverse representation p.l7|) in order to finally 
clarify the issue. 
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